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Abstract 

The mechanism of string creation for light quarks developed earlier is considered 
for nonzero quark chemical potential. A strong modification of the confining string 
due to finite quark density (chemical quark potential fi) is observed. As a surprising 
result in a multiquark system with a common string junction an attractive well 
appears of radius fj,/a and of an average depth equal to fi, which implies formation 
of multiquark hadrons. Preliminary estimates predict a new phase transition to 
multiquark hadron phase at rather high densities (in heavy ion collisions or in 
neutron stars) when neutron matter is compressed to 3-4 normal nuclear densities. 

1 Introduction 

Recently the high-density effects in QCD attracted a lot of attention because of the 
possible density phase transition [H [2] . On the fundamental level the physical expectation 
of any phase transition may be connected to the possible reconstruction of the vacuum and 
for that one needs that the energy density to be of the order of the vacuum energy density 
£cr ~ £vac ~ ~Y^c-^{{Fl^uy)y^cr ~ 1 GcV/fm'^ which provides the drastic change in 
the vacuum structure and may cause the phase transition. In case of temperature phase 
transition the corresponding energy density e is indeed of the order of Scr, and as it was 
first argued in [3] and later measured on the lattice [4], the QCD vacuum is strongly 
transformed in a way, that most part of colorelectric fields evaporate above T^. In this 
way Tc was calculated through e^r and found in good agreement with lattice data [5J, 
see [6] for a review. It was shown earlier in [7], that Chiral Symmetry Breaking (CSB) 
occurs due to confinement (more explicitly its colorelectric scalar components), which 
means that CSB should disappear simultaneously with confinement. In another way in 
the framework of hadron resonance gas and making use of the low-energy theorems at 
T 7^ a similar conclusion was obtained through application of the effective dilaton 
Lagrangian to gluodynamics and to QCD [8]. 

The same type of arguments for the density phase transition would imply that at the 
baryon density of ~ 1 nucleon/fm^, i.e. 3 -h 6 times higher than the standard nuclear 
density, the vacuum can be reconstructed in such a way, that part of fields, e.g. the 
colorelectric fields responsible for confinement, disappear above critical density. The ac- 
tual calculations in [5j have supported this conclusion and give critical chemical potential 
^g(T = 0) = 0.6 GeV. 
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This is the picture of deconfinement due to high density. However, of much more 
practical importance can be possible phase transitions in the confining region, at 2 ^ 3 
normal nuclear densities. 

It is the purpose of the present paper to start the investigation of the role of density 
on the vacuum fields in general and confinement in particular. 

In this paper we ask ourselves a short and simple question: how the nonzero baryonic 
chemical potential fi acts on the confinement of light quarks, and come to the unexpected 
answer, that the confining string of light quarks is destroyed gradually by /i in such a way, 
that the one part of string, near the string junction is eaten by the nonzero fi, while at 
distant r the string survives. We show that this can make multiquark states (MQS) more 
advantageous and thus cause a new type of phase transition. The plan of the paper is as 
follows. In section 2 the effective Lagrangian is derived from the QCD Lagrangian with 
nonzero n and their solution is discussed in section 3. Section 4 generalizes results to the 
case of baryons and multiquark bags. Section 5 is devoted to the physical implications of 
results and prospectives. 



2 Derivation of Effective Aq Lagrangian 

One starts with the QCD partition function in presence of quark chemical potential n in 
the Euclidean space-time, and we begin with the zero temperature, T = 0. 

where Sq{A) = ^ J(FI^^{x)yd'^x, m is the current quark mass (mass matrix rh in SU(3)), 
and the quark operator ^ilJaa{x) has flavor index a(/ = 1, ■■■Uf), color index a(a = 1, ...Nc) 
and Lorenz bispinor index a{a = 1,2,3,4), and we use the contour gauge [9j to express 
A^{x) in terms of F^^,. One has for the contour Zfj^^s, x) starting at point x and ending at 

Y = z{0,x) 

A,ix) = £ ds^^^^^^^F^Ms)^ = £dT,^,iz)F^,iz). (2) 
Integrating out the gluonic fields A^{x), one obtains 

Di/jDijj^e-f ^'/'"^(*^+*™-*^T4)-'^'^'='''^e^^^«^"^'^^«i"^'" (3) 
where the EQL proportional to {{A"^)) is denoted by L^eql^ 

L%L = jj d'xd'y /^Kx) /^(x) ^i,t,{y) ^^..(l/)(4ft^)4?(2/))7&^S^ (4) 
Average of gluonic fields can be computed using ([2]) as (see [TOj for details of derivation) 

/(4^(x)4^(i/)) = ^ r du,a,{u) I' dv,a,{v)D^'''''\u-v){5,Ji,-5,Ju,), (5) 

l\r Jo Jo 
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where D^^'^\x) is the correlator {Ei{x)Ei{0)) or {Hi{x)Hi{0)) . As it was argued in [TD] 
the dominant contribution at large distances from the static antiquark is given by the 
color-electric fields, therefore at the first stage we shall write down explicitly Lr^Qi^{el) for 
this case, i.e. taking = z/ = 4. As a result one has [TO] 

L%L{d) = ^Jd'xJ dS (x) (y) '^l^aeivW^^^^^ J\x, y) (6) 

where J^{x, y) is 

J^{x,y) = r du, r dviD^{u-v), i = 1,2,3. (7) 
Jo Jo 

One can form bilinears \E'^f = -^ipaa ^'^ae and project using Fierz procedure given isospin 
and Lorentz structures, ^'^f — > '^^^''^\x,y). Here we consider only ip^ij: bosonization. 
With the help of the standard bosonization trick (here J = -^J^) 

^ j ^^g^ jf/^DxeM-xh + ^'^h + ^x-m (8) 

Z = j D^D^+Dx exp Lqml (9) 
one obtains the effective Quark-Meson Lagrangian (QML) 

L^QML = jd'xj A { ^i^tai^Wd + tm- tii^,U5{x -y) + tM%'\x, y)] Va^Z/)- 

y)J\x, y)x^-^'\y, x)} (10) 

and the effective quark-mass operator is 

Mif{x,y) =Y.X^''''\x,y)Oi'^tfJj{x,y). (11) 

n,fc 

Here the operator O is a set of all irreducible combinations of Dirac matrixes. 

The QML in Eq. fllOl) L^qml contains functions x^"^'^^ which are integrated out in ([9]), 

and the standard way is to find x'^'^'^^ from the stationary point of Lq\.j^. Limiting oneself 
to the scalar and pseudoscalar fields and using the nonlinear parametrization one can write 
for the operator M in flTUl) 

M{x,y) = Ms{x,y)U{x,y),U = exp{i-f5(f)), ^{x,y) = (l)^{x,y)t^. (12) 

After integrating out the quark fields one obtains the ECL in the form 

4ci(M5,0) = -2nfN,{j''{x,y))-'Ml{x,y)+N,trlog[iid + tm-tfi^,)i + iMsU]. (13) 

^^{2) 

The stationary point equations —^ff^ = = at = 0O) = Mj°) immediately 

show that 00 = and M^^^ satisfies nonlinear equation 

tMl'\x,y) = AtrSJ''{x,y) = {^^S^,)j''{x,y), S{x,y) = -[id + tm - tfi^^ + iMsU]-l. 

(14) 

This equation plays the role of the gap equation and is the main point of our further 
investigation. For /i = this was done in [10] and in the next section we find how results 
of [To] are modified by the nonzero fi. 
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3 The confining string at nonzero /i 

Our basic equations ([6]), ([7]) are nonlocal in time because of the integral over dx^dy^ in 
(10). This nonlocality and the parameter which it governs can be handled most easily, 
when one uses instead of M{z, z'), S{z, z') the Fourier transforms. 

S{z, - z;, z, z') = / e^^'^(^^-<)^(p4, z, z')^ (15) 
and the same for M{z, z'). Then from (fTi|) one obtains a system of equations 

{pi — idz — im + z/i74)5'(p4, z,w) — i J M(p4, z, z')S'(p4, z', w)dz = 5^^\z — w) (16) 

To simplify matter, one assumes for D^{x) the Gaussian form, D^{x) = 



D{0) exp ( — ^ ) . Then for M(p4, z, w) one has 



^M(j94,z, w) = 2v/^T„ / ^e-(P*-P4)^7^lx 



X [J^(z,w)745(p'4,z,w)74] (17) 

where is defined in ([7]) and we have factored out the time-dependent exponent, using 
the Gaussian representation of D{u). 

All dependence of M onp^ as can be seen in (fT7|) is due to the factor exp[— (^4— P4)^T^] 
and disappears in the limit when Tg goes to zero, while the string tension a ~ D{0)Tg 
is kept fixed. This limit can be called the string limit of QCD, and we shall study its 
consequences for equations f|T6l) . f|T7j) in this section. 

So in the string limit, with M independent of p^, let us consider the Hermitian Hamil- 
tonian 

Hij^ = {j— + (3m- /x)V^„(z) +(3 j M{p, = 0, z, z')^„(z')rf^z' = eMM^) (18) 
with eigenf unctions ipn satisfying usual orthonormality condition 

llj:^{x)lp.m{x)d^X = 6nm, 



From f|T8l) it is clear, that one can redefine e„(/i) + fi = En, and and ipn do not depend 
on /i. Therefore in all subsequent formulas one can use the same equations as in [lOJ, but 
with the replacement e„ — e„ — /i. In particular, the Green's function S can be expressed 

as 

c/ N ^ 74^«(x)^+(y) . . 

'5p4,x,y =2^ 19 

n Pi- H^n - /i) 

Inserting (fT9|) into (fTTll one has integrals of the type: 

-^r^, ^ = T^liszgnie^ - /x)(l + 0{p,Tg, UT^) (20) 

-00 271 (p4 — i[en — /i)) 2 
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We are thus led to the following expression for M in the string limit 

M(p4 = 0, z, w) = v^Tg J^(z, w)74A(z, w) (21) 
where the definition is used 

^(z, w) = ^ i)n{'z)sign{en - /i)^/'^ (w) (22) 

n 

Let us disregard for the moment the possible appearance in M of the vector component 
(proportional to = 1,2,3,4) and concentrate on the scalar contribution only, since 
that is responsible for CSB and confinement. Then one can look for solutions of the Dirac 
equation (fTSjl in the following form [10] 



where /' = 2j — /, and introducing the parameter ^(j, /) = (j + ^)sign{l — j), and replacing 
M by a local operator (the generalization to the nonlocal case is straightforward but 
cumbersome, for a possible change in the nonlocal case see [ID]), we obtain a system of 
equations 

/ ^ + ^Gn-{en-fi + m + Mscai{r) - M,ect.{r))Fn = 

\ ^ - + (£„ - /i - m - Mscaiir) - M,ectXr))Gn = ^ ^ 

where we suggest M = M^^ai + liM-^ect. 

Eq.( l23i) possesses a symmetry (£„ — /i, F„, k) {ji — Sn, Fn, Gn, —n) which means 
that for any solution of the form ( l22l) corresponding to the eigenvalue — /x, there is 
another solution of the form 

corresponding to the eigenvalue (/i — £„). 

Therefore the difference, which enters (!22|) can computed in terms of Fn, Gn as follows 

A(z, w) = Ao(z,w) - AA(z,w) (26) 

where Aq is the value of A for /i = 0, i.e. the same as in [10], while AA is defined as 

AA(z,w)=2 J2 ^n(z)V^+(w). (27) 

0<£n</t 



Using decomposition (1251) one can write AA as 



A \(r- r'') — 9 \^ ( GnG^VtVf^, -iGnF^VlVL' 

ZAiHr,rj-z 2. \ ,FnGtn'Vt+, FnF+n'^+ 



0<e„,^i 



(2J 



where we have denoted Q = ^jim, ^' = ^ji'M, and we disregard nondiagonal part of A. 

At this point one can follow the relativistic WKB method for Dirac equation [1T\ 
applied to calculation of A in [10] in case of = 0. The classically available region for 
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ijjn{r) with energy e„ = e is r^in < r < r^ax, where rmax.min = , and the 

summation over n in ( 1271) transforms into integration over de, with the lower limit (for a 
given r) emin = or. In this way one has for the upper diagonal element in ( l28l) . 

2(j rfi/ar -t- _L 1 

AA(+,+) = — 5(1 -cos^rr') / ^ cos(aV7^^)g(/i - ar) (29) 

vr^r Ji — 1 

with a = crr|r — r'l, and we keep r ^ r' everywhere except for a, since for large a (when 
r is far from r') both A and AA fast decrease. 

In a similar way for the lower diagonal element in ( !28l) one has 



AA(-, -) = —5(1 - cos ^rr') / dr ^ cos{aV^^^)e{fi - ar) (30) 

TV^r Ji vr^ — 1 

and taking Aq in fl25p from [ID] the resulting form for A fl26p is 



Per , , f°° drcosi^ay 



A(r, r') = (3kscalar + "^Kector = ^"^(l - cos6'rr') / 



2cr _ „ . f/^/'^'" rdr 



-i^(5(l-cosV0 / cos(aVr2 - 1). (31) 

TT^r Ji A/r^ — 1 

Here /5 = 74, and 1 is the unit Dirac matrix, which means, that the second term on 
the r.h.s. of ( |3Tll contributes to the vector part of the resulting mass operator ( |2T1) . while 
the first term contributes to the scalar part. One should take into account, that in the 
first term, 

a ^.^ ^ . r°° (ir cos (a Vr^ - 1) ^ ^ 

Ascaiar = 1 " COS ^rr' / 7=^ 32 

TminifJ^) = fJ^/crr for > ar and 1 otherwise, so that for large r, r ^ ^, one has the 
standard //-independent value 

Kcaiar{r ~ r' > /i/cr) = -^Ko{a)S{l - cos6'rr') (33) 

where we have used relations for the McDonald function Kq 

, , cos axdx f°° -, , . , 

= / /T— / daK,{a) = -. 34) 
"'0 vl + Jo 2 

One can check that at large r,r'{r ~ r' > n/a) h^scaiar ~ ^icoZar ^ smeared 6 
-function, 

^ASi(r,r')rfV = l. (35) 



However for large /i, /i ^ AscaZar is different from AgcaiarifJ' = 0), and for r r' < fi/cr 
one has approximately 

A_,.,(r, r') = A^i - /(r, r')^(/. - ar) (36) 
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where 



f{r,r)= ; ^ , =/ ^===cosaA, (37) 



and Ao = ^(^)'-l. 

Let us now turn to the vector part of interactions, A^ector, 

2(J rfi/ar tcIt 

Kector = - cos6'rr') / ; ^ , cos(aVr2 - l)6'(/i - ar) = 

TT^r Ji Vt — 1 

9^ sin(aJf^)^ - 1) 
- i^S{l - cosOrr) ^^^^^ ^. (38) 

Returning back to Eq. ( !2TI) one can deduce, that 

M{p4 = 0, r, r') = v^Tg J^(r, r') [A^ca^r, r') + l4Avector] = M^cai + l4M^ect. (39) 

Taking into account, that at r, r' ^ Tg and for the Gaussian D^{x) one has from (JTj) 



rr' 



J^(r, r') = ^^2T.0FL)(O) min(r, r') (40) 



one has for M at r, r' ^ and for r = r' 

Mscai{r,r') = (Tr(p)(r,r') -e(r,rO) (41) 

where 6^^\v,r') = A^(!^°-'(r, r') and 

e(r,r') = 5^ = - cos9rr')f{r,r')9{fi - ar) (42) 

and f{r,r') is defined in fl37|) . 
For M^ect one has, using ( l38i) . 

M,eet(r, r') = -2ar^^{r, r')9{fi - ar) (43) 

where 

^ sin(aA/f^)^ - 1) 

(^^(r,r') = — 5(l-cosV0 ^^A^^:^ . (44) 

7r"^r a 

Note, that one should actually symmetrize all these expressions, e.g. ^ etc., 

but we always are in the regime, where r ~ r'. To estimate the magnitude of nonlocal 

kernel M5caz(r, r') and M^,erf(r, r') it is convenient to introduce in (|T8|) the local limit of 

the mass operator, namely 

Mscal,vect{r) = J dh' M seal ,vect{r , T ) . (45) 

Exploiting the equalities 

J dh'd^^'^r, r') = I dh'^ir, r') = J dh'^^{r, v') = 1 (46) 

one arrives at the expressions 

Mscai{r) = are{ar - ii), M^ect = -2ar9{fi - ar). (47) 

In the next section we shall discuss approximations made in deriving pij) . fjl3l) . P7|) 
and physical implications of these results. 
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4 Multiquark states and baryons at nonzero density 



Using the local limit for interactions of a light quark with string junction, Eq. (H71) . 
one can calculate the change of masses of white states due to density. At this point it is 
necessary to define exactly the physical meaning of the parameter /i, introduced in (ITU]) as 
the quark chemical potential. That would be true, in the deconfined phase, when quarks 
can move freely over the volume and occupy the Fermi sphere in momentum space. In 
the confined phase, however, the term fi is meant to be attached to the common string 
junction of all quark states, constituing the complete set of states ( !24|) . This meaning is 
especially clear from Eq. (l2Ti) . where fj, exemplifies the boundary of the occupied states 
{e} in the {n,j, n} notation. 

Therefore, assuming that multiquark states located at different string junction points 
do not overlap, one can associate /i for a given quark with the largest value of occupied 
by other quarks, belonging to the same string junction point. In what follows we shall 
show, that in some situations it will be advantageous for two or more nucleons (with dif- 
ferent string junctions) to coalesce into a common state with one string junction, building 
in this way the contracted potential (jTTj) . 

Note, that for one nucleon (3 quarks) the contraction mechanism does not work and fi 
can be taken at zero value, since no forbidden states exist for each quark of a given color. 
This, again, is true, provided one-quark states do not overlap at the given density. The 
situation changes, however, for two nucleons, since e.g. one can add to the first nucleon, 
i.e. 3 quarks in S-states, 3 quarks in P-states relative to the same string junction point. 
The latter will move in the potential (jUj), where /i is equal to the S-state Dirac eigenvalue 
Es, and the energy ep{fi = Es) may become less than es{fi = 0): 



The same type of inequality may occur for higher L states, when more nucleons coalesce 
to the same Multiquark State (MQS). 

This would imply instability of nuclei with respect to transition to MQS matter. 

If, however, Eq. (1481) does not hold in the limit of zero density, one may have an 
interesting situation for increasing density, large baryonic chemical potential fiB- Note, 
that fiB 7^ 3/i, since /i is attributed to MQS string junction and does not grow with fiB 
unless different MQS start to overlap. Hence for > fj48l) should be replaced by 



To estimate, let's take es{fi = 0) ^ 500 MeV, ep{fi = 500 MeV) ^ 630 MeV (see the 
Tables 1,2 below), so for (Fermi) ?a 130 MeV the phase transition can occur (which 
corresponds to densities 3 -h 4 times higher than a nuclear one). 

Below we perform calculations of e{fi) for different values of /i and L = 0, 1 We start 
with the values of en,L,K{0), which were computed earlier numerically and analytically 
using WKB approximation [10] (see Table 1). 

In Table 2 we show the change of EniK as AM5 = 3Ae„(S') for S and AMp{3q) = 
2Aen{S) + Aen{P) for the P-wave baryon due to nonzero value of /i, AsnLnifJ') = ^nLK^Ij) — 
^nLniS^), calculated in perturbation theory, considering change of potential in (j47j) as a 



ep(/i = Es) < ssifi = 0) 



(48) 



SpifJ' = £s) <es{li = 0)+ e_B (Fermi). 



(49) 
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Table 1: Eigenvalues ^ computed numerically, as in fl2^ 



1 





1 


2 


n 


siguK 









1.29 


1.93 


2.40 





+ 




2.02 


2.51 


1 




2.34 


2.76 


3.11 


1 


+ 




2.79 


3.18 



Table 2: Mass shifts of 3g system due to nonzero Hq 



liq (MeV) 


AM5(3g) MeV 


AMp(3g) MeV 


100 


-5 


-4 


200 


-70 


-50 


300 


-266 


-200 


400 


-605 


-483 


500 


-1023 


-876 



perturbative parameter. One should note at this point, that the vector part of interaction 
written in (H7j) for the case of light quark with heavy antiquark, and for the case of 3 
quarks interacting with string junction the vector part is transformed as an addition to 
the color Coulomb potential with coefficient 1/2 (see [11] for explicit derivation of scalar 
and vector interaction in the nucleon). This prescription was used for calculations in 
Table 2. 

5 Discussion of results 

Results of the previous section Eqs. fHTj) . fH3l) . fHTl) . can be formulated as follows. The 
relativistic WKB analysis leads to the /^-dependent modification of the confining string, 
where the piece [0,/i/(T] of the string near the origin of the string (situated at the heavy 
quark position in case of heavy-light quark, or at the string junction position in the case 
of baryons), is dissolved, and the linear confinement starts beyond the critical radius 
^cr = /i/o". Moreover, an attractive vector interaction appears in the same interval with 
the average magnitude (My^ct) ~ Z^- 

As shown in Table 2 the change in quark energies Ssifi), £^p(/u) is large and Aep{fi = 
400 MeV ) = -80 MeV, Aep{fx = 500 MeV) = -190 MeV, while \ep{fx = 0) - ^^(/i = 
0)1 ~ 320 MeV. Hence it seems to be not advantageous for a quark in neighboring nucleon 
to add in the P-wave to the 3q system in S state, with a reduced confinement. 

However, the gap is not large. Indeed, the MQS with five quarks in S state and one 
quark in P state will have the mass shift A only 130 MeV above the NN mass. Note, that 
phenomenological value of the 6q MQS A = 210 MeV is in the same ballpark Hence 
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Figure 1: Interactions Mscai{r) (solid line) and M^ectir) (dashed line)as functions of dis- 
tance r, with b = ix/a. Dotted line and dashed line show the qualitative smoothed form 
of both terms respectively. 
V 




for higher density this value of 130 MeV can be compensated by baryon Fermi energy, 
and the baryon matter can become unstable for transition into MQS matter. 

These conclusions should be taken as qualitative. First of all, the WKB method is 
not a good approximation at small distances, and we have omitted exponentially damped 
part of V'n(z) in the spectrum, therefore the inner part of the string is to some extent 
delocalized (see [lO] for details) and smoothed, as shown in Fig.l. 

Secondly, the interaction between MQS and baryons was not considered above. 

Thirdly, we have not taken into account a possible modification and destruction of 
the vacuum due to the influence of high density quark matter, which might decrease a 
or cancel the string completely (as it is happens in the thermal phase transition [6]). 
We have calculated in [5] the critical value of /ig(T = 0), where deconfinement due to 

density occurs /i^"* = ^'^^'^^ + (^^^2) = 0.6 GeV, which corresponds to 7 times 
normal nuclear density. We are interested here however in the confinement region and 
small temperatures, with density few times larger than normal. 

If however, no density induced vacuum deconstruction takes place, then the resulting 
physical picture according to Eqs. (H71) . is the net decreasing of confinement in the inner 
region of some ensemble of quarks, and appearance of attractive vector potential of the 
order of /i acting on each quark. This may cause creation of partly deconfined bubbles 
consisting of 3n quarks, n = 2, 3, .. in the midst of the nuclear medium, and dynamically 
is similar to the 3ng bag formation, which was studied before in the framework of the 
Quark Compound Bag model [H]. Note, however, that bag boundary conditions might 
be strongly modified as compared to the standard MIT bag model. As it is, preliminary 
estimates in the previous section demonstrate a possible formation of a new phase of 
MQS at rather high densities (in heavy ion collisions or in neutron stars) when neutron 
matter density exceeds 3-4 normal values. A quantitative analysis of this situation needs 
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a more accurate consideration of the 3n quark system using nonlinear equations for the 
3nq Green's function, generahzing Eq. f|T^ . 

The formation of these high-density 3nq bubbles may be connected with the explana- 
tion so-called cumulative effects in the hadron-nucleus(and nucleus- nucleus) collisions, [15J 
for an example of this discussion see [I6j and refs. therein. 
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